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Abstract The invention of the laser immediately en- 
abled the detection of nonlinear photon-matter interac- 
tions, as manifested for example by Franken et al.'s de- 
tection of second-harmonic generation. With the recent 
advancement in high-power, high-energy lasers and the 
examples of nonlinearity studies of the laser-matter in- 
teraction by virtue of properly arranging lasers and de- 
tectors, we envision the possibility of probing nonlinear- 
ities of the photon interaction in vacuum over substan- 
tial space-time scales, compared to the microscopic scale 
provided by high-energy accelerators. Specifically, we in- 
troduce the photon-photon interaction in a quasi-parallel 
colliding system and the detection of higher harmonics 
in that system. The method proposed should realize a 
far greater sensitivity of probing possible low-mass and 
weakly coupling fields that have been postulated. With 
the availability of a large number of coherent photons, 
we suggest a scheme for the detection of higher har- 
monics via the averaged resonant production and decay 
of these postulated fields within the uncertainty of the 
center-of-mass energy between incoming laser photons. 
The method carves out a substantial swath of new ex- 
perimental parameter regimes on the coupling of these 
fields to photons, under appropriate laser technologies, 
even weaker than that of gravity in the mass range well 
below 1 eV. 



1 Introduction 

Recent astronomical observations suggest the existence 
of non-luminous stuff in the universe: dark matter and 
dark energy. Understanding their origin is one of the 
greatest scientific challenges of the 21st century. It may 
be natural to expect that even the universe is just a 
collection of vacua around us. Therefore, we consider 
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a novel method to directly probe the structure of vac- 
uum under laboratory conditions that relate to these 
issues. To our current knowledge no direct signature of 
the non-luminous stuff has ever been observed in ter- 
restrial experiments. This is in spite of the fact that 
numerous advanced theories exist, including axions [I], 
minicharged particles [2,i3ji4n5j, and dark energy [6j[7]. 
The reason for this may be thought of as follows. After 
Rutherford's discovery of the inner core (i.e. nucleus) of 
an atom being very tiny compared to the dimension of 
the already tiny size of the atom, the experimental search 
went to explore ever smaller constituents of matter and 
thus the thrust went for higher-energy or momentum 
experiments. Theories have gone hand-in-hand with this 
exploration, succeeding in ever shorter-ranged interac- 
tion theories and unification of forces, as exemplified by 
the electroweak theory [8j . We refer to this standard and 
extremely successful method as the high-momentum ap- 
proach. Almost all laboratory research efforts have been 
on this approach to date. Although successful in explor- 
ing high energy physics, this approach is not suitable 
to explore energies much lower than 1 eV. These fields, 
that might exist in much lower domains than 1 eV, can- 
not strongly couple to matter, because if they did, they 
would have been long ago discovered in the low-energy 
region. Thus these fields, if they ever exist, must couple 
weakly. This means that we need an extremely strong 
driver to manifest a sufficiently strong signal overcom- 
ing this weakly coupling interaction, showing up above 
noise. So far such sufficiently powerful photon sources 
did not exist. However, this may be changing now with 
more intense lasers becoming available |^ . What we have 
called the high-amplitude or high-field method jTUlfTT] 
may provide an alternative path to detect such low- 
mass, weakly interacting fields that are spread over semi- 
macroscopic scales. 

When the energy of the constituent is much lower 
than that of the probing photon, conceptually we may 
employ two laser beams in a co-propagating geometry. 
The two coUinear beams produce a very low center-of- 



K. Homma", D. Habs and T. Tajima 



mass energy interaction at the beating frequency (being 
equal to the difference of the two laser frequencies) [12] . 
This interaction could resonate with the very low eigen- 
frequency of the constituent, should there be an eigen- 
mode in its vicinity. Second of all, the co-propagating 
setup allows us to make the two beams interact over a 
much prolonged interaction time, thus much amplifying 
the nonlinearities and signals arising from these. 

In order to pick up the experimental signal of a 
strong coupling to the long-range mode, we suggest 
using higher-harmonic generation. The pioneering re- 
search by Franken et al. 13i detected the nonlinear- 
ity in a quartz crystal via second-harmonic generation. 
Two photons in co-propagation accentuate the interac- 
tion through the quartz fields over the coherence vol- 
ume in order to initiate second-harmonic generation. 
This process may be schematically looked upon as the 
case displayed in Fig. [ija). There the quartz nonlinear- 
ities mix two forward-propagating photons (uj) to pro- 
duce a photon with 2a; (and possibly another photon 
with frequency ^ 0, also referred to as optical rectifica- 
tion or difference frequency mixing). The quantum elec- 
trodynamics (QED) process is illustrated by Fig. db). 
Two incoming photons are mediated by virtual electron- 
positron fields and outgoing are two photons. The ex- 
treme forward-scattering amplitude with quasi-parallel 
incident photons is known to be largely suppressed in 
the QED process, as we discuss later. This is because 
the center-of-mass energy of the colliding two photons is 
too low to satisfy the relevant mass scale of the electron- 
positron pair. 

Given these hints from analogous pictures on how to 
probe medium-like features of vacuum, this paper rather 
considers a simple photon-photon scattering process via 
the averaged resonant production and decay of light- 
mass fields within the energy uncertainty in the center 
of mass system between incoming laser photons. This 
gives a solid basis for the design of experiments based 
on intense co-propagating lasers, as much independent 
as possible of models on the light-mass fields in vacuum. 
As we discuss later in detail, the process we focus on is 
based on Fig. [Ijc). 

The QED process as illustrated in Fig. \V[h) is de- 
scribed by the Euler-Heisenberg effective Lagrangian il4i 
in the low-frequency limit of two incoming photons. This 
predicts that the vacuum under the influence of an elec- 
tromagnetic fleld induces a birefringence characteristic. 
The ratio between the first and the second term in the 
Lagrangian can yield a general test to see whether the 
vacuum contains other effects beyond QED, examining 
its value at 4:7. In general, a scalar field (p and a pseu- 
doscalar field a in vacuum may contribute to the first 
and second term, respectively. Light scalar fields as can- 
didates of dark energy have been recently intensively 
discussed [5], while the pseudoscalar fields (axion-like- 
particles) may be a source of dark matter fP and also 
possibly dark energy [TB] . We suggested how to measure 
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Fig. 1 Schematic diagrams of photon-photon interactions 
in matter and in vacuum, (a) Second harmonic generation 
in the experiment by Franken et al. [13., arising from the 
nonlinearity of crystal fields irradiated by (intense enough) 
laser fields; (b) Probing QED vacuum nonlinearities as sug- 
gested by Heisenberg and Euler |14II15] . where the vertices of 
the coupling in the Feynman diagram are characterized by 
the fine structure constant of the vacuum a. Thus a weak 
nonlinearity requires much more intense fields than in the 
case (a). The leading order interaction is the elastic photon- 
photon scattering, though as a higher order there exists a 
second harmonic generation as well; (c) Probing potential 
light-mass m fields in vacuum with intense laser fields. The 
vertices are characterized by very feeble couplings of 
and g [7|- The expected second-harmonic generation may be 
said to be not different to case (a). In order to increase the 
observable signal, suggestions have been made. 



the phase shift of a probe laser across an intense elec- 
tromagnetic field, based on the phase-contrast Fourier 
imaging [T7]. However, we recognize that the test by the 
phase-contrast Fourier imaging is rather limited in the 
mass and coupling of those fields [TH] . Therefore, we ex- 
tend our method to search for those new types of fields by 
instituting co-propagating laser beams. This approach 
may be looked upon as in Fig. [IJc). Again two paral- 
lel photons come in, while two parallel photons come 
out. As we discuss later, the frequency shift in vacuum 
is simply explained by the strong Lorentz boost in the 
quasi-parallel colliding system, which causes the blue- 
and red-shift of photon energies emitted to the forward 
and backward directions, respectively. We note that our 
approach is similar to, but distinct from many laboratory 
experiments with lasers p^[20ll2n[22l[23ll24l[25] already 
performed and proposed to search for those low-mass 
fields. 

Consider more details of the effective interaction La- 
grangian L as illustrated in the triangle part in Fig.[IJc), 
where the scalar field (j) and the pseudoscalar field a cou- 
ple to the electromagnetic field via quantum anomaly- 
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type couplings [7] 



(1) 
(2) 



where F^i^ — dA^/dx" — dA^/dx^ is the antisymmet- 
ric field strength tensor with the four-vector potential 
Af^ of the electromagnetic wave and its dual tensor 
F^"' = l/2e^"'"'F,^, with the Levi-Civita symbol e^'"^". 
The difference (f> versus cr is in their allowed couplings 
to polarization states of two photons as we discuss later. 
Here gM~^ provides the coupling strength, g and M 
carry the subscripts (fi and cr, respectively, indicating 
the corresponding type of fields. The dimensionless cou- 
pling g is typically proportional to the dimensionless fine 
structure constant a for the two photons to couple to the 
virtual charged-particle pair in the triangle part. The ef- 
fective coupling includes the large mass scale M to cou- 
ple to the light field with a mass m. The large M induces 
the weakness of the coupling via M^^. For example, the 
Newtonian constant G is expressed as SttG = hcMp^, 
where Mp is the Planckian mass of 10^^ eV. The weak- 
ness of G is the manifestation of the large mass scale at 
the vertex in the triangle coupling. In what follows we 
omit the subscripts (j) and a on the coupling gM~^ and 
the mass of the light field m, unless we need to explic- 
itly distinguish the type of the fields. We use the natural 
units h = c = 1 throughout the subsequent sections, 
unless explicitly noted. 

As a quite challenging case we have attempted a the- 
oretical approach to search for an extremely light scalar 
field as a candidate of dark energy via the averaged reso- 
nance scattering process in [53] . Given intense laser fields 
in the near future, the method may provide a new win- 
dow into scoping physics on the Planckian mass scale 
by photon interactions in a quasi-parallel incident laser 
beam in the laboratory. In this paper we review the 
essence of the approach and further develop basic for- 
mulae in Sect. [5] to apply the method to a specific exper- 
imental setup considered in Sect. [31 in order to discuss 
reachable limits on the coupling strength gM~^ and the 
mass m for a laser intensity attainable within the ELI 
project [17] in Sect. [3| In Sect. [5] we provide a perspec- 
tive on how our method serves as a new type of scope 
to probe the semi-macroscopic vacuum, which has not 
been extensively investigated yet. 



2 Quasi-parallel photon-photon interaction in 
low-mass and weakly coupled fields 



central line along the z axis. We define the z ^ x plane 
formed by pi and p2. The components of the 4- momenta 
of the photons are given by pi = (w sinz?, 0, a; cos t9; 
and the same for p2 but with the sign of ?? reversed, 
and p3 — (wa sin ^3, 0, a;3 cos 6*3; ws) and p4 with 0^3, 6*3 
replaced by 104,— 64, respectively. The angles 9^ and 64 
are defined as shown in Fig. [2] This coordinate system 
can be transformed to the center-of-mass (CM) system 
for the head-on collision (z9 = 7r/2) by the Lorentz trans- 
formation with v/c —5- 1 for i9 — > 0. Conversely, this im- 
plies that the realization of the quasi-parallel collision 
in the laboratory frame corresponds to the realization of 
an extremely low CM energy, as we see below. 

In this frame one of the final photons in the for- 
ward direction along the z axis must have an upshifted 
frequency due to the energy-momentum conservation, 
independent of the physical origin of the dynamics. In 
the limit of — ?> 0, a process of UJ3 2uj is realized. 
This frequency doubling nature is an extremely valu- 
able characteristics from the experimental point of view, 
as compared to the case with no frequency shift in the 
center-of-mass system. In addition, more importantly, it 
is essential to maintain a quasi-parallel nature of the in- 
cident beams in order to access resonance, as we shall 
stress later. 

The energy-momentum conservation laws requires 
following relations; 

0-axis : UJ3 + uj^ = 2uj, (3) 
z-axis : CJ3 cos ^3 + uj^ cos 9^ = 2uj cos i?, (4) 
X-axis: 0^3 sin (?3 = 0^4 sin ^4 . (5) 

From the conditions < ^3^4 < 2a;, we may choose 
O<03<?9<04<7r, without loss of generality. From 
Eq. ©-([SJ we derive the relation 



sm 6^3 = sm 6^4 



sin^ {} 



1 — 2 cos d cos 04 -I- cos^ -d 



(6) 



The differential elastic scattering cross section per 
solid angle dfi'i favoring the higher photon energy uj^ is 
given by 



^^{8T:uj)-^sin-^^L03/2iof\M\^, (7) 

ttJ ^3 



where Ai is the invariant amplitude and 



2.1 Kinematics oj quasi-parallel system 

As illustrated in Fig. [2l we introduce an unconventional 
coordinate frame, in which two photons labeled by 1 and 
2 sharing the same frequency are incident nearly paral- 
lel to each other, making an angle d with the common 



UJ SIT? 1} 



LJ3 



1 — COS -d cos 6*3 



(8) 



Here we expect the upshifted frequency W3 — >■ 2uj, as 
03 for — >■ 0, as mentioned before. 



4 



K. Homma", D. Habs and T. Tajima 



Pz e(i) 




Fig. 2 Definitions of kinematical variables for tlie suggested 
co-propagating piiotons (tliis figure is quoted from [26]). 



2.2 Dynamics of two-photon interaction via resonance 

We consider the scattering amplitudes only for the case 
when light-mass fields are exchanged via resonance (s- 
channel amplitude) . The resonance decay rate of the low- 
mass field with the mass m into two photons is expressed 
as 



m 



(9) 



The low-mass field is exchanged between the pairs 
and (p3,P4), thus giving the squared four- 
momentum of the field 



9' = (Pi +P2) = (cos2i9- 1) 



(10) 



with the metric convention (+ + H — ) for the definition 
of the four momenta. Qs corresponds to the CMS energy 
of the photon-photon collision. 

With the polarization vectors given by ej'^\ where 
i = l,---,4 are the photon labels, whereas (3 ~ 1,2 
are for the kind of linear polarization as depicted in 
Fig. [5J we summarize the non-zero invariant amplitudes 
for scalar field exchanges 

A^iiii = A^2222 = -A^ii22 = -X2211, (11) 
and for pseudoscalar field exchanges 

A^1212 = A^1221 = -A^2112 = -A^2121, (12) 

where the first two digits in the subscripts correspond 
to the states of the linear polarization of the incoming 
two photons 1 and 2, respectively, and the last two cor- 
respond to those of the outgoing two photons 3 and 4, 
respectively, as illustrated in Fig. [2] 

We focus on one of these non-zero amplitudes by de- 
noting it a,s A4; 



uj^ (cos2i?- 1)^ 
2w2(cos2?9- 1)-Hm2' 



(13) 



where the denominator, denoted by V in the following, 
is the low-mass field propagator. We note that in (|T0)) 
is time-like. We then make the replacement 



(14) 



W and 



(m — ir) w — 2imr. 

Substituting this into the denominator in Eq. 
expanding around m, we obtain 



2? « -2 (1 - cos 2i?) ix + ia) , with x = - w^, (15) 



where 



mV2 



l-cos2i?' l-cos2i? 
From Eq. ^ and (ITCl) . a is also expressed as 



(16) 



(17) 



which explicitly shows the proportionality to Af"^. We 
then finally obtain the expression for the squared ampli- 
tude as 



(18) 
is largely 



As for the off-resonance case x ^ a, 
suppressed due to the factor a? cx M"* for the case of a 
small coupling M^^ . On the other hand, if experiments 
take the limit of a; ^ Wr, \M.\'^ — >■ (47r)2 is realized 
from (fT8)) ideally. This is independent of the smallness 
of the factor M~^, as expected from the off-resonance 
case or equivalently from the square of p^ . This is the 
most important feature arising from the resonance that 
overcomes the weak coupling stemming from the large 
relevant mass scale such as M = Mp. However, we are 
then confronted with an extremely narrow width a for 
e.g. < 1 eV , M ~ Mp = eV and 1 eV. 

We now discuss how to overcome this difficulty. 



2.3 Enhancement by averaging including resonance 

In conventional high-energy collisions, the beam mo- 
menta are implicitly supposed to be their mean values. 
This is because the momentum spread, or the uncer- 
tainty expected from the de Broglie wavelength of the 
relativistic particle, is negligibly small compared to the 
relevant momentum exchanges in the interaction that ex- 
periments are interested in. Resonance searches in such 
experiments must adjust x in (|15D and Eq. (ITS)) 
such that the mean x is close enough to the peak lo- 
cation within ±a. On the other hand, in the case of co- 
propagating laser beams aiming at the detection of ex- 
tremely small momentum exchanges via the resonance, 
the situation is quite different due to the nature of the 
incident waves. This is because the uncertainty included 
in the initial photon momenta is much larger than the 
relevant energy scale of the resonance, leading to the 
condition |x| ^ a. In this case the squared scattering 
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amplitude must be integrated over the possible uncer- 
tainties on the incident wave function on an event-by- 
event basis. As we will explain below, the uncertainty 
in the incident photon momenta is related to the uncer- 
tainty in X via the uncertainty in the incident angle 'd. 
Therefore, it is instructive to consider the feature of the 
integral of the resonance function of the Brcit-Wigner 
(BW) formula [28] from x- to x+ as follows; 

^= / atan^ , (19) 

where / = a7r/2 and I — air for x+ = ~X- = and 
X+ = —X- = OO; respectively. This indicates that the 
value of the integral is proportional to a, i.e., M^^ from 
Eq. (|17p . The value ranges for the finite and infinite inte- 
grals over only a factor of two. From this fact, we expect 
that the integral enhances the squared scattering ampli- 
tude by a factor of compared to the non-resonant in- 
teraction proportional to M~'^ from Eq. (fT5|) . as long as 
the peak is contained within the experimental resolution 
on X, *-e., the condition x-i- > a and X- < is satis- 
fied. This implies that experiments in the co-propagating 
laser beam configuration need no efforts to adjust x close 
to the extremely narrow resonance region, thanks to the 
subsequent huge enhancement by the integral over the 
wide range on x- Meanwhile, it is difficult to identify the 
exact location of the resonance mass within the wide gate 
on X- Since we are interested in having a sensitivity to an 
extremely weak coupling such as gravity M = Mp, the 
enhancement of the squared amplitude is more crucial 
than finding the exact location of the resonance masses. 
As we discuss in the following sections, however, we may 
be able to provide a crude estimate on the order of the 
mass scale of the resonance even in such a situation. 

The consideration above leads us to parametrize the 
squared scattering amplitude as follows: 

/ + 00 
p{pi,P2)\M\^dx 
-OO 

= {Anf / p{p,,p^) r—^dx, (20) 

where p{pi,p2) is the normalized probability distribu- 
tion to supply nominal combinations (j5i,P2) from the 
incident laser fields. The real part x{pitP2) is indirectly 
specified by {pi,P2) via the incident angle d between the 
incident two photons. This parametrization expresses av- 
eraging over possible combinations of pi and p2 . We note 
that {pi,P2) are not a priori the observed momenta, but 
just a nominal specification among the possible initial 
momenta. In other words, p{pi,p2) is not a statistical 
weight on the discrete momenta after the contraction of 
the wavepacket of each photon state. This implies that 
an infinite statistics is not necessary to obtain the con- 
tinuous nature. Rather, we need this treatment even for 
a two-photon state as long as the source of photons is 
not a perfect plane wave. This allows for a continuous 



integral on x via the continuous combination of p(pi,p2) 
in Eq. ((20|) . As long as the probability weight p{pi,p2) 
is close enough to unity around the resonance peak, the 
enhancement discussed with Eq. (|19l) is guaranteed. This 
is the essence of our main strategy of the co-propagating 
configuration in order to overcome the difficulty due to 
the narrow resonance width a. 

In order to design experiments, we start from the 
resonance condition, the first of Eq. (jl6p . by assuming 
< 1, 

m ~ 2l}r^^r (21) 

where the subscript r in both angle and energy refers to 
a state satisfying a resonance condition. We note that 
the product 2'drUJr, corresponds to the CM energy of the 
incident two photons. This indicates that experiments 
have two adjustable handles for a given mass scale or 
the CM energy. We emphasize that we can lower the 
CM energy by several orders of magnitude by only in- 
troducing smaller values of i? with a fixed w. This should 
be contrasted to high-energy colliders, where a large ef- 
fort is needed to increase the CM energy by an order 
of magnitude. This advantage is also supported from a 
technical point of view, since scanning the incident angle 
d should be much easier than scanning the energy uj of 
the resonance. We point out that the resonance condi- 
tion in Eq. (j2ip is not just given at one point, but rather 
in a hyperbolic band in the iS — lj plane, with a finite 
resolution with 5'd of the incident angle 'd. This implies 
that the deviation Suj from the resonance energy can 
satisfy the same resonance condition with a different 
within ±S-d. As far as Suj/ojr ^ Sd/dr is satisfied in the 
setup, we can ignore the effect of 5uj. This is in fact the 
case, as can be seen in the following discussions. There- 
fore, we can take the attitude that we fix the incident 
energy at the optical frequency u — ujopt and scan m 
by changing i!) around tJ^, where ujopt and i^r satisfy the 
resonance condition based on Eq. ((2T|) 

<t = mV(4t?^). (22) 

From Eq. (PT|) and Eq. we see that detecting lower- 
mass fields requires smaller values of 'dr, namely the par- 
allelism of incident photons. In this case the variation on 
leads to the variation on x- This is expressed by the 
following relation based on the second part of Eq. ((TS]) 

X(t?) = <t - ^ = - (23) 

where e = •d/'dr in the unit of is introduced. 

We now discuss the average of the squared amplitude 
|A^|2 over the possible uncertainty on the incident angle 
I? 

|7W|2 / p(d)\M\^d^ 




p{s)\M\^drde, (24) 
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where p{'&) is a probability distribution function nor- 
malized between and n/2 as a function of the continu- 
ous uncertainty on z9 from arbitrarily chosen two-photon 
combinations within a laser pulse. The incident angle 
•d is re-expressed with e by dd — drde for the second 
part of Eq. (gH). From the relation Eq. we find 

e = (1 - x)-^/"^ and de = l/{2Lolp^)e'^dx. Equation ^ 
is then re-expressed with % 



(1 - xf'^ ix' + 



dx, (25) 



where Eq. ([18)) is substituted. This equation is the exact 
representation of Eq. (|20p , starting from the uncertainty 
on the incident angle "d, if the upper limit of the integral 
range is regarded as large enough compared to a. Let us 
define a; = to explicitly discuss the structure of the 
integral kernel in units of the width a of BW. With ^, 
Eq. (pSj) is further re-expressed as 

1?, 



a-Hl-(2^.A) = } p((l_a^)-l/2) 



(1 - a^yi^ e + 1 



dL (26) 



where the first factor of the integral kernel corresponds 
to a normalized weight function and the second is BW 
with a width of unity. This expression explicitly shows 
the enhancement by the factor of a, implying the pro- 
portionality to M^^ based on Eq. (IT71) . As long as p is a 
monotonic function, the weight function in front of BW 
can be close to unity for small ^ because of <C 1. With 
such a weight we expect that the value of the integral 
may be close to that of BW, as we discuss with a specific 
weight function in the following section. 

The remaining issue is how to further cope in ex- 
periments with the problem of still very small values of 
Af~^, although much larger than M^^. First, this can 
be solved by the sin"'* d behavior of the cross section in 
Eq. (O, that arises from the phase volume factor and the 
flux factor in the quasi-parallel two-photon interaction. 
For an extremely light mass, this factor gains a large 
number due to the small -dr- Second, the intense laser 
fields can provide a large luminosity and the intensity of 
the signal is proportional to the square of the intensity 
of the laser in the limited case of an incoherent two- 
photon interaction. We have three ingredients or knobs: 
the enhancement by the weighted BW integral, the 
"dr dependence and the growth of the laser intensity. By 
marshalling these knobs, we expect to increase the de- 
tectability for undiscovered low-mass fields in vacuum, 
which have evaded from our grasp to date. 

In the following sections we consider experimental 
realizations with a;^ ~ 1 eV (optical laser), aiming at 
the mass range as low as possible. We then plug explicit 



weight functions into Eq. based on the suggested 

experimental setup. By combining Eq. (1261) and Eq. ([7]), 
we discuss reachable mass-coupling limits for a given 
laser intensity attainable in near future experiments. 



3 Second-harmonic detection in the 
Quasi-Parallel System 

We emphasized the importance of the photon-photon 
interaction in a quasi-parallel system or a small-angle 
setup, in order to enhance the signal due to low-mass 
constituents. A simple way is to explore the mass range 
m < TTw by using two independent laser beams with a 
small incident angle. We then directly measure the res- 
onance curve in Eq. (I2ip by scanning both and lo to 
quantitatively observe the nature of the resonance curve. 
For the much smaller mass scale, or equivalently smaller 
incident angle, however, we must take into account the 
beam spread in the diffraction limit. This determines the 
controllable smallest incident angle, or the mass scales of 
the light fields, we look for. We consider here the case of 
a single focused laser beam, in order to provide the sim- 
plest basis to quantify reachable mass-coupling limits for 
a given set of experimental parameters. We concentrate 
on the detection near the second-harmonic, rather than 
the laser frequency itself, to enhance the sensitivity of 
detecting the photon-photon interaction. 

The conceptual experimental setup with a single- 
beam focusing geometry is illustrated in Fig. [3] Incident 
photons from a Gaussian laser pulse with linear polariza- 
tion are focused by a ideal lens into the diffraction limit. 
Quasi-parallel incident photons interact with each other 
between the lens and the focal point, from which photons 
3 and 4 are emitted in nearly opposite directions along 
the z axis with ~ 2a; and a;4 ~ 0. The dichroic mir- 
ror is transparent for the non-interacting photons with 
the beam energy of w, while W3 is reflected to the prism 
(equivalent to a group of dichroic mirrors), which selects 
W3 among residual w and sends it to the photon detector 
placed off the 2;-axis. This process is assisted by a polar- 
ization filter. From the polarization dependence of the 
invariant amplitude in Eq. (|TT|) and Ec^. (fT2)) . the com- 
binations of polarizations of two photons between the 
initial and final states must satisfy 11 —7> 11(22) for a 
scalar field exchange and 12 — > 12(21) for a pseudoscalar 
field exchange, respectively. We note that we can choose 
the type of fields we search for by setting the initial po- 
larization state. In the case of single-beam focusing, the 
search for a scalar field is easier, because we do not have 
to mix the two polarization states as in the case of a pseu- 
doscalar field. Furthermore, the selection of the rotated 
final state 22 can enhance the signal-to-background ratio 
for the scalar field case, because a huge number of non- 
interacting photons has the final polarization state of 
11. In what follows we provide formulae to evaluate the 
accessible limit on the mass-coupling defined in Eq. ^ 
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for a given laser intensity, in the case that we detect a 
double- frequency photon per laser shot. 

The Gaussian profile is a basic constraint in typi- 
cal laser fields, where the aperture of a lasing material 
has a finite size in the transverse area. The solution of 
the electromagnetic field propagation in vacuum with 
a Gaussian profile in the transverse plane with respect 
to the propagation direction z is well-known |29j . The 
electric field component in spatial coordinates (x, y, z) is 
expressed as 



E{x,y,z) oc 

Wq 



w{z) 



exp < ~i[kz — H{z)] — r 



wizf 2R{z) 



where k — 27r/A, r — \/ + J/^, is the minimum 
waist, which cannot be smaller than A due to the diffrac- 
tion limit, and other definitions are as follows: 



w{z) 



Wo 



1 



R^z 1 



H{z) 



■y2 



tan"" ( — 

.ZR 



ZR 



X 



(28) 



(29) 



(30) 



(31) 



Based on the Gaussian laser parameters above, we now 
estimate the effective luminosity C over the propagation 
volume of the laser pulse. We now restore the physical 
dimensions of h and c in this section, unless explicitly 
noted. Consider a Gaussian laser pulse with duration 
time r, with the speed of light c and an average number 
of photons N per pulse. The exchange of a low-mass 
field may take place anywhere within the volume defined 
by the transverse area of the Gaussian laser times the 
focal length / before reaching the focal point. We first 
consider the effective number of photons Nint during an 
interaction with the time scale of As a result of the 
interaction we observe a frequency-doubled photon in 
the laboratory frame. The momentum transfer of ~ huj/c 
between photons defines the minimum interaction time 
scale from the uncertainty principle as follows 



At > 2ttuj" 



(32) 



The effective number of photons during At is expressed 
as 



At - 

N,nt = —N. 
T 



(33) 



Making the pulse duration r ~ At maximizes the instan- 
taneous luminosity. Suppose a point z along the laser 



R = oo 



P2 



?1\ 



P1 

F 




polarization ©jj 
filter ^ 



photon detectorO 



Fig. 3 Suggested experimental setup for the co-propagating 
photon interaction and detection. The linear polarizations of 
incident and outgoing photons are drawn only for the scalar 
exchange with a scattering amphtude IAI1122I as an example. 



propagation axis. The instantaneous luminosity at the 
point z is defined as 



C{z) 



C{N„ 



t:w'^{z) 



2TrwQ -|- zj^ 



(34) 



where C{Nint, 2) denotes a combinatorics to choose two 
photons amongst a large number of photons available 
within the time scale At, and the expression w^(z) in 
Eq. is substituted to obtain the second part with the 
approximation for the combinatorics. We then consider 
the averaged instantaneous luminosity C over the focal 
length / as follows. 



/ 



C{z)dz • 



~-tan-i(//zfl,) 



27r/u;^ 



zr^ 



Ar2 

"*tan-i(//zfl) 



2/A 



(35) 



The number of effective bunches h is related with / as 



zAt' 



(36) 



The effective luminosity C over the propagation volume 
of the laser pulse is finally expressed as 



f Nf 



cAt 2/A 



tan-i(//zfl) 



At 

~2ctX 



tan-i(//z^)(37) 



where Eq. ([55)) is substituted in the last step. 

The minimum beam waist wq s-t z = varies with 
the experimental conditions, the focal length /, and the 
diameter of the incident beam d. 



Wo 



nd ' 



(38) 
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where the approximation is vahd when /A ^ {Tr/A)(P. 
The uncertainty on the incident angle between two hght 
waves is then expected to be 



A 

TTWo 



d 



(39) 



We see that Ad is controlled via wq by choosing suitable 
values for / and d in experiments. 

The possible uncertainty on the incident angle § 
affects the average of the squared amplitude \M\'^ as 
shown in the first part of Eq. In order to obtain 

an approximation close enough to reality, we plug the 
following step function into Eq. (IMl) : 



l/Ad forO<z9<Z\?? 
for Zi?? < I? < 7r/2 



(40) 



which is normalized to the physically possible range < 
■d < tt/2. By substituting Eq. (gO]) into Eq. we 
obtain 



\M 



2 — 



(47r)^ dr 



a-^{l-{-»^/A^f} 



(1 - aO^/' (C' + 1) 



(41) 



where the first factor of the integral kernel corresponds 
to the weight function and the second is the Breit- 
Wigner function (BW) with a width of unity. The weight 
function of the kernel is close to unity for small ^, 
due to the smallness of a in Eq. (flT)) . Therefore, the 
value of the integral in Eq. (j4ip is almost equivalent to 
that of BW [30]. This is because the monotonic posi- 
tive weight function approaches zero as ^ — oo more 
rapidly than the pure BW, whereas the pure BW sup- 
presses the increase of the weight function close to zero 
at C ^ a-i{l - {'&r/A'&)^} for Ai} < 1. We then ap- 
proximate Eq. (|4T|) as the integrated BW over ± oo 
as follows: 



(47r)^ §r 



(42) 



Let us introduce the integrated form of the cross section 
in Eq. ([7]) over the solid angle 27r sin O^dO^ up to 9^ , which 
is equivalent to the lower limit of the frequency uJ^/lo = 
2 — h though the relation of Eq. ^ by specifying the 
frequency deviation < h ^ 1 from the exact 2uj as 
follows 



aih) 



\M\^ 



inujy sin i9 Jo ^ ^ 



(—] 2TTsm63d93 



TT{h/2){8TTUj)-^d-^\M\^, (43) 

where i? <i; 1 is used for the approximation in the last 
step. With a in Eq. ^ and JM^ in Eq. (gl]), the 



weighted cross section in Eq. is finally expressed 



"""-sIt) [it, )(w ) <«' 

where the approximation z?^ ^ 1 is also taken into ac- 
count. 

Multiplying Eq. 1^ by Eq. (gl]), we obtain the in- 
tegrated yield y above the lowest frequency 0^3 specified 
by h per laser pulse focusing as follows: 



y = Ca{h) 
1 At A 
256 T cT 



Ai} J \ M I 



There are several experimental knobs to affect the ob- 
servable events in Eq. (HSI) . If we choose r ~ Z\i ^ A/c, 
resulting in ct ^ A, we can maximize the effective lumi- 
nosity. From Eq. (1391) , the reduction of A-d or increasing 
wq enhances the yield in the case of very low-mass par- 
ticle exchange. From Eq. (1311) and (l38|) . we express //z_r 
as 

//«~.^. m 

From this relation, a shorter focal length enlarges 
f / zji. This introduces a slight increase for the factor 
tan~^(//zfl), though its effect is tempered by the na- 
ture of tan^^. 

As a short summary, we make the most important 
note from the experimental point of view based on this 
conceptual design. The condition dj./ A^} = 1 maximizes 
the chance to search for a resonance, while dr/A-d > 1 
results in a huge suppression of the cross section by 
{Ti = c = 1) as we discussed. This is because the reso- 
nance peak is out of the region covered by A-d. This 
parameter corresponds to a sharp cut-off of the cross 
section. Therefore, controlling A-d via the relations of 
Eq. dsn), Eq. ([35]), and Eq. can provide an ex- 

perimental way to define the mass range that we ex- 
clude, if no signal is found. From the resonance condi- 
tion in Eq. (|21l) . we evaluate the cut-off on the mass 
range we can exclude by this conceptual design as fol- 
lows {h = c = 1): 



nicut = 2A-duj. 



opt 



f 



to. 



opt I 



(47) 



where Eq. ([M]) was substituted. We note that the mea- 
surement is still sensitive to the mass range even below 
this cut-off, however, we cannot state where the mass 
is located below this limit. Therefore, if we cannot find 
any symptoms below this cut-off, we can exclude the en- 
tire mass range below this cut-off. On the other hand, 
if something is found, we must make effort to lower the 
cut-off by introducing a longer focal length and a smaller 
beam diameter in order to determine the location of the 
mass at which the signal disappears. 
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4 Sensitivity enhancement by high-intensity 
lasers 

As a demonstration we now discuss how much high- 
intensity lasers improve the accessibility to the weak cou- 
pling domain for the following reference case, based on 
a quasi-parallel colliding system. 

First, we note that At in Eq. (1321) is the resolvable 
minimum time scale. As long as we discuss an extremely 
low-mass field, the interaction time scale may be over 
h/mc^. In the case of low mass below 1 eV, the inter- 
action time scale may be much longer than 27ra;^(. On 
the other hand, photon-photon interactions must occur 
within a laser pulse duration time r. This fact implies 
that we should assume At/r = 1 even if At ^ r. Since 
the beam waist of focused laser becomes smaller as it 
approaches the diffraction limit, a short-pulse laser with 
the duration time close to maximizes the the ef- 

fective luminosity per laser shot. Therefore, we assume 
this condition in the following discussion. 

We are now ready to estimate the sensitivity to the 
coupling g / M for a given average number of photons N 
in the laser pulse, in the case that y is found to be zero 
per laser shot based on Eq. (|45p for the basic experimen- 
tal parameters discussed above: ujopt ^ 1 eV, r ~ 10 fs, 
At/r = 1, d ^ 1 m, and /i = 0.1 to require the frequency 
close to 2ujopt- Figure m shows accessible lower bounds 
on g/M as a function of m via the search for higher- 
harmonic generation by focusing a single laser shot. The 
blue dotted and solid lines indicate cases for N = 2 J 
and TV = 20 kJ, respectively. The upper and lower lines 
on each line style correspond to the case for / = 1 m and 
/ = 1 km, respectively. As we discussed with Ea. (|T7)) . 
the cut-off values on mass are lowered by introducing 
longer focal lengths. The sky-blue band shows the resolv- 
able mass range by gradually changing focal lengths from 
1 m to 1 km in which we can determine a mass based on 
the disappearance of higher harmonic signals from the 
appearance state. The black shaded area is the excluded 
region by conventional laser-driven experiments so called 
"Light Shining through a Wall (LSW)" [S]- In LSW a 
laser pulse propagates under a static magnetic field B 
and a low-mass particle a is assumed to be produced via 
two photon coupling to the particle (cr — i3— laser cou- 
pling). Since a is expected to couple with matter only 
weakly, a can penetrate a wall which prevents penetra- 
tions of laser photons on the other hand. The particle 
may couple again to i?-field placed over the wall result- 
ing in a photon in the induced decay process. 

Our arguments so far have been based on the ap- 
proach in which photons both in initial and final states 
are treated incoherently, giving an observable yield pro- 
portional to 7V^|7Wp. It is worth noting, however, that 
the intrinsic nature of laser fields may further improve 
the sensitivity to even smaller couplings as weak as 
the gravitational coupling g/Mp ^ 10~^''GeV~^ with 
g a = 1/137 which is indicated by the magenta line 



in Fig|31 Degenerated photons in a laser beam can in- 
duce the decay of resonantly produced low-mass fields 
into a degenerated final state as in the case of the static 
magnetic field to induce the decay of low-mass parti- 
cle used for LSW. If this is the case, compared to their 
spontaneous decays in vacuum, we can expect an en- 

1/2 

hancement factor ^ N^^^ in the scattering amplitude, 
caused by the creation operator to the degenerated state, 
where Nind is the number of degenerated photons in- 
volved in the induced laser field [5^126) . For example, 
if we limit h in Ea. (j43p to a narrow range around 0.5 
in the spontaneous decay process, a phase space decay- 
ing into CJ3 l.duJopt and LU4 ^ O.bujopt can be selected. 
If we induce the UI4 ~ O.buiopt by the degenerated state, 
the probability to emit W3 ~ 1.5a; opt photons is enhanced 
accordingly. If we rescale all photon frequencies by a fac- 
tor of two, what is suggested here is that we simply look 
for higher-harmonic generation via the following process: 
2cj + 2a; loj + 3a;, where la; is assigned as the induced 
field. In this case the yield of the 3a; wave should fol- 
low N'^ Nind\-M\'^ ■ Similar ideas based on the classical 
treatment are discussed in [551151] . The red lines in FigH] 
show the lower bounds when N = Ni„d is assumed. The 
conditions and conventions are same as those for blue 
lines except this additional induced field. It is a bit of 
a surprise that it is not impossible to reach weaker cou- 
pling domains even beyond the gravitational coupling in 
a scattering experiment by accumulating statistics over 
a reasonable time period, if we could introduce the in- 
duced luj waves into the same geometry as 2uj to make 
them focus at a time. We note that this method is also 
applicable to nsec pulses, typical for MJ-class laser fa- 
cilities, as long as the energy per pulse is larger than 
that of short laser pulses which may compensate the 
longer time duration by the N'^Nind dependence of the 
higher harmonic yield. The actual coupling limit must 
be evaluated, eventually based on the statistics of the 
background of higher harmonics. 

A major instrumental background for the frequency 
shifted radiation is in principle expected to be higher- 
harmonic generation (HHG) from the final focusing op- 
tical element and the one to reflect HHG to the photon 
detectors. The dominant source of HHG may be the in- 
terface between the residual gas and the surface of the 
optical element, where the centrosymmetry is maximally 
broken. Even from the maximal estimate ^ lO^^W/cm^ 
for a typical damage threshold which is a much lower 
intensity compared to HHG due to relativistic motion 
of surface electrons, we expect a negligible amount of 
10^^" for the case of second-harmonic photons from a 
1 m^ aperture size with a lOfs irradiation, if the op- 
tical components are housed in a vacuum containing 
10^° atoms/cm^ 10"^ Pa) [5S]. The conflrmation of 
a negligible contribution of the background from higher 
harmonics is a crucial subject for the present concept. 

As a dominant physical background we expect the 
lowest-order QED photon-photon scattering with a for- 
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10"^^10""*10"^^10"'^10-"10"'° 10"^ 10"° 10-^ 10"^ 10"^ 10"^ 10"° 10-2 10-^ 1 

m reVl 



Fig. 4 Accessible mass-coupling domain by searching for higher-harmonic generation. Blue and red lines are cases for 
spontaneous and induced decays, respectively. The dotted and solid lines indicate cases for the mean number of photons per 
laser pulse N — 2 J and iV = 20 kj, respectively. The upper and lower lines on each line style correspond to the case for 
focal length / = 1 m and f — I km, respectively. The sky-blue band shows the resolvable mass range. The shaded area is the 
excluded domain by laboratory laser-driven experiments "Light Shining through a Wall (LSW)"[3T]. See text for the detail of 
the other experimental parameters. 



ward cross section ~ (a^/TO^j^w^z?^ [35]. This turns out 
to be much smaller than Eq. (H^ . due to the specific 
behavior with respect to the incident angle i). This indi- 
cates that the lowest-order QED contribution is negligi- 
ble. 

If there will be no signal in the single-beam focusing, 
we only have to update the condition such that it satisfies 
A-d > -dr for heavier masses by increasing A-d. In such 
a heavier mass region, however, the two-beam crossing 
geometry relaxes the constraints on the optical design 
such as the focal length. In either case the single-beam 
focusing setup considered in this paper provides a basis 
to define the mass-coupling limit as well as the necessary 
beam intensity as we have demonstrated here. 



5 Conclusion 

We have suggested an approach to probe the nature of 
vacuum with intense lasers. The resonance search for 
extremely light-mass fields via higher-harmonic genera- 
tion has been explored by focusing high-intensity laser 



beams. This is similar to the idea already developed to 
probe matter. 

In this method we take note of the nonlinearities of 
vacuum that are speculated to wait for our sensitive de- 
tection of an extraordinarily feeble signal. In order to 
detect weak nonlinearities, we need to spectacularly en- 
hance the signal. The large leap in enhancing these sig- 
nals is garnered by the combination of (i) the rapid devel- 
opment of the intense laser technology and its adoption 
here; (ii) the employment of our suggested scheme allows 
for an enhanced interaction with the pursued fields. The 
former element (i) may be brought in, for example, by 
an intense optical laser beyond 1 kJ. For the latter factor 
(ii) we have suggested a method to include a resonance 
with co-propagating photons for the exploration of pos- 
sible new low-mass fields, with the aid of induced laser 
fields to promote the decay into a degenerate vacuum 
state. 

With the detection of higher-harmonic generation in 
the co-propagating setup, we should be able to survey a 
large sweep of the energy domain of the intermediating 
vacuum fields. If and when we pick up some signal in 
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Fig. 5 Experimental domains of various approaches to probe matter and vacuum as a function of the system size vs. the 
energy density. Selected systems are LC: electron-positron collision in the center-of-mass energy Ecms ~ 1 TeV at the future 
linear collider [57] i assuming the electron size 10~^*cm which is the upper limit of the electron radius obtained in high-energy 
collider experiments; LHC: proton-proton collision in Ecm = 14 TeV at the Large Hadron Collider [37]; RHIC: gold-gold 
collision in Ecm = 200 GeV per nucleon pair at the Relativistic Heavy Ion Collider [37], the rest proton indicated by the 
asterisk as the origin of this plot; ELI: an optical laser pulse expected in the ELI project [27] ; Schwinger: the Schwinger 
limit [39], Sun, the Milky Way Galaxy and the cosmic horizon with Otot ~ 1-0 and h = Ho/WO [km/s/Mpc]~ 0.7 [38]. The 
energy density axis is qualitatively interpreted as the inverse of the force range or the mass scale m of the exchanged force, 
because the mean free path becomes shorter in higher density states, as long as the coupling to matter is not weak. On the 
other hand, the coupling strength to matter gM~^ as defined in Eq. ([T} qualitatively reflects the necessary size of matter or 
vacuum in order to make the interaction manifest. The arrow to the higher energy density towards the Schwinger limit is the 
direction to probe nonlinear QED interactions and also towards an understanding of the non-perturbative nature of the intense 
field. The arrows directing to the lower energy density region indicate the extensible domain by using co-propagating intense 
laser fields, since the sensitive mass range is below ~ 1 eV covering mass scale relevant for dark energy and the coupling may 
be probed to a scale as weak as gravitational coupling for lighter mass scales. The energy density in this direction depends 
on the context. In the context of the scalar field as a candidate of dark energy in [6ll26]. the energy density should be close to 
that of the cosmic horizon. 



one particular energy range, perhaps we can zoom in to 
this specific energy (and thus wavelength) of photons by 
arranging the various knobs, such as the crossing angle 
and the (long) beating wavelength of the electromag- 
netic waves. By setting up a specific resonance cavity, 
we may be able to further increase the sensitivity and 
more deeply study their properties. 



Given a highly intense optical laser beyond 1 kJ per 
fs-pulse duration in the near future, the realization of 
these suggestions may become an exciting challenge for 
future experiments investigating the physics of the vac- 
uum. We might be able to reach coupling strengths as 
weak as the gravitational strength by utilizing induced 
laser fields. 
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Figure E] illustrates the experimental domains of var- 
ious approaches to probe matter and vacuum, in terms 
of the system size as a function of the energy density. 
The energy density axis is qualitatively interpreted as 
the inverse of the force range or the mass scale m of 
the exchanged force in Eq. ([T^ . because the mean free 
path becomes shorter in higher-density states as long as 
the coupling to matter is relatively strong. On the other 
hand, the coupling to matter gM~^ in Eq. (IT5|) or Eq. ([T]) 
qualitatively reflects the necessary size of matter or vac- 
uum in order to make the interaction visible. 
The Galileo-type telescope observes gravitational phe- 
nomena. These are on the extremely weak coupling scale 
of Mp^ with zero mass exchange. High-energy particle 
colliders, the Rutherford-type microscopes, focus on par- 
ticle generation phenomena. These are due to strong cou- 
plings with heavy mass exchanges within the fm scale. 
There is a huge gap between these two approaches. In 
other words, the region of weak couplings with finite but 
light mass exchanges has hardly been probed so far. It is 
quite natural to start exploring if there exist important 
pieces of the puzzle of nature in these domains. These 
explorations might grant us deeper understanding of the 
nature of vacuum such as dark energy [40j . The progress 
of modern physics has been simply driven by those two 
types of experimental approaches. The proposed method 
with high-intensity lasers probes the semi-macroscopic 
vacuum compared to particle physics and on a much 
smaller scale of vacuum compared to cosmology. Pro- 
vided such semi-macroscopic vacuum scope, we increase 
our observational window into a new parameter regime 
of the vacuum. 
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